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Abstract 

We study biased, diffusive transport of Brownian particles through narrow, spatially periodic 
structures in which the motion is constrained in lateral directions. The problem is analyzed under 
the perspective of the Fick-Jacobs equation which accounts for the effect of the lateral confinement 
by introducing an entropic barrier in a one dimensional diffusion. The validity of this approxima- 
tion, being based on the assumption of an instantaneous equilibration of the particle distribution 
in the cross-section of the structure, is analyzed by comparing the different time scales that char- 
acterize the problem. A validity criterion is established in terms of the shape of the structure and 
of the applied force. It is analytically corroborated and verified by numerical simulations that the 
critical value of the force up to which this description holds true scales as the square of the peri- 
odicity of the structure. The criterion can be visualized by means of a diagram representing the 
regions where the Fick-Jacobs description becomes inaccurate in terms of the scaled force versus 
the periodicity of the structure. 

PACS numbers: 05.60.Cd, 05.40.Jc, 02.50.Ey 
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I. INTRODUCTION 



In many transport phenomena such as those taking place in biological cells, ion channels, 
nano-porous materials and microfluidic devices etched with grooves and chambers, Brownian 
particles, instead of diffusing freely in the host hquid phase, undergo a constrained motion. 
The uneven shape of these structures regulates the transport of particles yielding important 
effects exhibiting peculiar properties. The results have implications in processes such as 
catalysis, osmosis and particle separation y, 0, 0, Q, 8, [q] and for the noise- induced 
Tansport in periodic potential landscapes that lack reflection symmetry (ratchet systems) 
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12j or ratchet transport mechanisms which are based on asymmetric geometries. 
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termed "entropic" ratchet devices |10l. Illl 
separation of DNA fragments in narrow channels 
shape. The translocation of structured po 
to determine their sequence and structure 



13|. For example, it has been found that the 
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16j is largely influenced by their 



ynucleotides through nanopores also allows one 



17|,ll8|,ll9(. 



The motion of the particles through these quasi-one- dimensional structures can in prin- 
ciple be analyzed by means of the standard protocol of solving the Smoluchowski equation 
with the appropriate boundary conditions imposed. Whereas this method has been very 
successful when the boundaries of the system possess a regular shape, the challenge to solve 
the boundary value problem in the case of uneven boundaries represents typically a very 
difficult task. A way to circumvent this difficulty consists in coarsening the description by 
reducing the dimensionality of the system, keeping only the main direction of transport, but 
taking into account the irregular nature of these boundaries by means of an entropic poten- 
tial. The resulting kinetic equation for the probability distribution, the Fick- Jacobs (F-J) 
equation, is similar in form to the Smoluchowski equation, but now contains an entropic 
term. The entropic nature of this term leads to a genuine dynamics which is very different 



from that observed when the potential has an energetic origin [20|. It has been shown that 
the F-J equation can provide a very accurate description of entropic transport in 2D and 
3D channels of varying cross-section 20, l21|. 

However, the derivation of the F-J equation entails a tacit approximation: The particle 
distribution in the transverse direction is assumed to equilibrate much faster than in the 
main (unconstrained) direction of transport. This equilibration justifies the coarsening of 
the description leading in turn to a simplification of the dynamics, but raises the question 
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about its validity when an external force is applied. To establish the validity criterion of a 
F-J description for such biased diffusion in confined media is, due to the ubiquity of this 
situation, a subject of primary importance. 

Our objective with this work is to investigate in greater detail the F-J approximation for 
biased diffusion and to set up a corresponding criterion describing its regime of validity. We 
will analyze the biased movement of Brownian particles in 2D and 3D periodic channels of 
varying cross section and formulate different criteria for the validity of such a F-J description. 

The paper is organized in the following way: In Sec. [Tll we describe the physical situation 
and introduce the model defined through the corresponding Langevin and Fokker-Planck 
equations. In Sec. Illlt we introduce the F-J approach for the unbiased situation and extend 
it to the driven case. Sec. [IV] is devoted to establish a criteria for the validity of the F-J 
approximation derived by comparing the different characteristic time scales of the process. 
In Sec. lYl the accuracy of the F-J description is tested against numerical simulations for a 
2D periodic channel, and the conditions of validity are summarized in a diagram in terms 
of the scaled force versus the periodicity of the structure. In Sec. I VII we provide further 
explanations on when and why the equilibration assumption fails and the F-J approach leads 
to useable results. Finally, in Sec. IVIII we present our main conclusions. 

II. DIFFUSION IN CONFINED SYSTEMS 

In typical transport processes through pores or channels (like the one depicted in Fig. 
[1]), the motion of the suspended particles is induced by application of an external potential 
V{r) resulting in a force F. In general, the dynamics of the suspended Brownian particles 
is governed by Langevin 's equation: 



where r*is the two or three dimensional position vector of a particle of mass m, rj is its friction 
coefficient, k-B the Boltzmann constant, T the temperature, and a dot over the quantity refers 
to time derivative. The thermal fluctuations due to the coupling of the particle with the 
environment are modeled by a zero- mean Gaussian white noise ^(t), obeying the fluctuation- 
dissipation relation {Ciit) ^jit')) = 2 6ij6(t — t') for z,j = x,y,z. In the over-damped case. 
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FIG. 1: Schematic diagram of the channel confining the motion of the biased Brownian particles. 
The half-width w is a periodic function of x with periodicity L. 

i.e. when mr{t) « rif{t), the inertia term in Eq. ([T]) can safely neglected and the Langevin 
equation describing the dynamics of a Brownian particle within the channel reads: 

^ ^ = -VV^(f(t)) + Vvk^m ■ (2) 

In addition to Eq. ([2]), the full problem is set up by imposing reflecting boundary condi- 
tions at the channel walls. 

The corresponding Fokker- Plank equation for the time evolution of the probability dis- 



tribution P{f,t) takes the form 
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23j: 



dt 

where J{f,t) is the probabihty current: 



^^(-*)^-V/(,<). (3) 



/(r, i) = - + £>oV I P(f, i) , (4) 



and 



Do = kBT/ri , (5) 

denotes the diffusion coefficient of the suspended particles. 

Due to the impenetrability of the channel walls, the normal component of the probability 
current J{f, t) vanishes at the boundaries. If n denotes the vector perpendicular to the 
channel walls, the reflecting boundary conditions read: 



J(f, t)-n = 



(6) 



In this paper we focus on the case of a symmetric 2D channel where the force is constant 
and directed along the axis, cf. Fig. [H The half-width of the 2D channel is given by a 
periodic function lj{x), i.e. uj{x + L) = uj{x) for all x. In this case, the boundary condition 
reads 



dco'(x) 
dx 



F ^ dP{x,y,t) 
-P{x, y, t) - Do 

7] OX 



dP{x,y,t) _ 

J^O 7\ — U 

oy 
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aX y = -^uj{x). For an arbitrary form of uj{x), the boundary value problem defined through 
Eqs. ([3]), (jlj) and (JTj) is very difficult to solve. Despite the inherent complexity of this 
problem an approximate solution can be found by introducing an effective one-dimensional 
description where geometric constraints and bottlenecks are considered as entropic barriers 
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III. THE FICK-JACOBS APPROXIMATION 



In the absence of an external force, i.e. when F = 0, it was shown 2J, |25| that the 
dynamics of particles in confined structures (such as that of Fig. [1]) can be described by the 
F-J equation 



dP{x,t) 
di 



d f d P 

9x \ ° dx h{x) 



(8) 



obtained from the 3D (or 2D) Smoluchowski equation after elimination of y and z coordi- 
nates. Here P{x, t) is the probability distribution function along the axis of the 2D or 3D 
channel, h{x) is the dimensionless transverse cross-section h{x) := tt [uj{x)/L]'^ in 3D, and 
the dimensionless width h{x) := 2uj{x)/L in 2D, where uj{x) is the radius of the channel 
in 3D (or the half-width of the channel in 2D). This description is valid for |ci^'(a;)| <^ 1 
(the prime refers to the first derivative, i.e. uj'{x) = duj{x)/dx), and it has been shown that 
the introduction of an effective x-dependent diffusion coefficient can considerably improve 
;he accuracy of the kinetic equation, thus extending its validity to more winding structures 
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261]. The expression: 



D{x) 



(1 + cj'(a;)2)" 



(9) 



where a = 1/3, 1/2 for two and three dimensions, respectively, has been shown to accurately 



account for the curvature effects 



2l|. 



In the presence of a constant force F along the direction of the channel the F-J equation 
can be recast into the expression [20|, [21 1 



dt dx \ ^^dx ^ k^T dx J ^ ^ 

which defines the free energy A{x) := E — T S = —Fx — T k-Q ln/i(x), where E = V = 
—Fx denotes the energy contribution and S = ks lnh{x) the entropy contribution. For a 
symmetric channel with periodicity L, the free energy assumes the form of a periodic tilted 
potential. 

Note that Eq. ffTOj) may typically describe the time-evolution of a ID particle distribution 
within an energy landscape. In the present context, however, due to the reduction of the 
geometric confinements in 2D or 3D space into one dimension, we end up with an entropic 
contribution to the free energy. In absence of a fixed bias F = 0, we deal with a pure 
entropic situation with the free energy reading A{x) = -Tk-B lnh{x). Equation ( ITOi) then 
corresponds to a F-J equation with a spatial-dependent diffusion coefficient. Likewise, for a 
situation involving solely transport in a one- dimensional energy landscape, the free energy 
reduces to A{x) = V{x) with duj{x)/dx = 0, yielding to the well-known Fokker-Planck 
equation in ID, with a constant diffusion coefficient Dq, reading 

Recently, we have shown that the dynamics of a confined Brownian particle, in presence of 
an applied bias, can accurately be described by means of equation ( ITOl) 20l |. In the presence 
of very strong applied bias, and for more winding structures, however, the F-J equation 
becomes inaccurate. In the present work we present further numerical and analytical results 
and will set up tailored criteria under which the F-J approximation assumes good validity. 



A. Particle current 



One of the key quantities in transport through quasi-one-dimensional structures is the 
study of the average particle current (x), or equivalently the nonlinear mobility, which is 
defined as the ratio of the average particle current and the applied force F. For the average 
particle current we derive an expression which is similar to the Stratonovich formula for 
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the current in titled periodic energy landscapes [27|, l28|, l29|, l30| , but with a spatial diffusion 
coefficient. A detailed derivation of this expression is given in the Appendix|Xl cf. Eq. flA12l) . 
Hence, we obtain the nonlinear mobility for a 2D channel with a shape defined by u{x): 

1 l-exp(-/5FL) 



1 _|_ f dujjx') ^ ^ 
^ \ dx' 



[12] 



where we have made use of Eq. (Q and the relations Dq = k^T/i] and (3 = l/k^T. Substi- 
tuting z' = x'/L and z = x/L, the nonlinear mobility scaled with the friction coefficient r) 
can be expressed, for the case of constant forcing, in terms of a single, dimensionless scaling 
parameter {2^: 



f:=PFL. 

Therefore, Eq. ( fT2i) for A{x) = —Fx — (3~^ \n{h{x)) leads to 

/.N 1 l-exp(-/) 
/ Jo d2/(z,/) 



(13) 
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where 



I{z,f)=exp{fz)h{Lz)J^ dz'exp{-fz')h~\Lz') 1 + 



lduj{Lz'] 
L dz' 



(14b) 



Eq. (IT^ can be transformed into Eq. (6) in Ref. 20(] by interchanging the order of integration. 

The asymptotic values of the nonlinear mobility for the cases / — *^ and f ^ 00, can be 
evaluated to read: 

1 



lim /u(/) 1] = 

{h{x)){^h-^ix) 



(15) 



and 



lim /i(/) r] 

f^oo 



Do 
D{x) 



(16) 



where 
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g{x)dx , 



(17) 



is the average over a spatial period, given an arbitrary periodic function g{x). 



IV. VALIDITY OF THE FICK-JACOBS DESCRIPTION IN PRESENCE OF A 
CONSTANT BIAS 



The reduction of dimensionality done implicitly in the formulation of the F-J equation 
relies on the assumption of equilibration in the transverse direction. An estimate of the 
conditions under which equilibration occurs can be made by analyzing the different time 
scales involved in the problem. For the sake of simplicity, let us focus on the situation of a 
2D channel, although the same discussion can readily be extended to 3D. In a 2D channel in 
the presence of an external force, in the axial direction, one can identify three characteristic 
processes, with the corresponding different time scales. One is diffusion in the transverse 
direction over a distance A?/, whose time scale is 

Ay' 



(18) 



Similarly, the time scale associated to diffusion in the axial direction is 



The third time scale is defined through the characteristic time associated to the drift 
(ballistic motion) over a distance Ax given by 

= (20) 

Tdrift V L 

where we have used the scaling factor / = FL/k-BT. Rearranging the previous expression, 
we obtain 

LAx , . 

Tdrift = ■ (21) 

In order to have a good equilibration in the transverse direction, the characteristic time 
scale associated to diffusion in this direction has to be much smaller than the other two 
time scales. Therefore, in the absence of an external force, equilibration in the transverse 
direction occurs if Ty/T,j. <^ 1. This results in the condition: 

(22) 

which constitutes the validity criterion of the F-J approach such as put forward by Zwanzig 
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In presence of a force along the axis, equilibration in the transverse direction demands 
that the condition — ^ <^ 1 also holds. Consequently, 

J y ^ J ^ > < 1 ^ (23) 



2LAx L2 

where in the second step we have replaced the characteristic distances Ay by the width 
2c(j(x), and Ax by L. 

A general estimate of the criteria that has to be satisfied is that max(— , — ^) ^ 1. An 

^ ^Tx^ Thrift ^ 

even stronger criteria in order for the F-J description to hold in presence of a constant force 
can be put forward by considering the sum of the two ratios, cf. Eqs. ( !22|) and (!23|) . i.e.: 

Eq. (I24p provides a quite stringent criteria that indicates when the F-J description of 
a system is expected to be valid. Note also that this is a local criterion, i.e. for a given 
channel, there would be regions (associated to drastic changes in the shape of the channel, 
i.e uj'{xY ^ 1 ) where equilibration in the transverse direction is not feasible, whereas in 
others is fulfilled. It is then more convenient to work with a global criterion of validity rather 
than with a local one. One way of getting that global condition is by averaging the local 
criterion over the period L of the channel, yielding one of our main results: 

(c.'(x)2) + |^(c.(x)^)«l. (25) 

In order to get an explicit estimate of the dependence of the maximum force value on the 
periodicity L of the channel, we define a critical force value /c, for which the inequality fl2S]) 
becomes an equality, i.e. {uj\xY) + jj{uj{xY) = 1. Then the critical force value reads: 

fc = 7^7^, (1 - i^'i^r)) , (26) 

2{uj{xj'') 

which indicates that the critical force scales asymptotically as L^, if we fix the overall shape 
of the channel and change only its periodicity. 

Eq. ( |26l) provides an estimate of the minimum forcing beyond which the F-J description 
is expected to fail in providing an accurate description of the dynamics of system. The 
quantitative value of the critical force will obviously depend on the level of accuracy sought. 
What is really important is how it depends (or scales) with the relevant parameters of the 
problem. 
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FIG. 2: (color online) The numerically simulated (symbols) and analytically calculated, cf. 
Eq. (jl4p . (lines) dependence of the scaled nonlinear mobility /u(/)?? vs. the scaling parameter 
/ = FL/k^T is depicted for two channels in 2D with different spatial periods. For both chan- 
nels the scaled half width is given by uj{x) = sm{2'irx / L) + 1.02; L = 1 : diamonds and solid 
line (blue), L = 27r : circles and dotted line (green). The dashed line indicates the deterministic 
limit ^{f)ri = {x)/{F/r]) = 1. 

V. NUMERICAL SIMULATIONS FOR A 2D CHANNEL 

In order to check the consistency of the criteria proposed and the validity of the F-J 
description in the presence of a force, we will compare the analytical result for the scaled 
nonlinear mobility obtained in Eq. (IT^ with the simulation results of the overdamped 
Langevin dynamics in Eq. ([2]) for the 2D periodic channel sketched in Fig. [H We remark 
here that the extension of our scheme to 3D with a rotational symmetry along the transport 
axis is possible as well. This will consume more computation time, but the overall findings 
remain qualitatively robust. This feature has been verified with a few numerical tests. 

The shape of the channel is described by 



where a controls the slope of the channel walls and 2(6 — a) gives the channel width at 
the bottlenecks. In order to use dimensionless quantities we refer any physical length to 
a, i.e. the scaled boundary condition then reads uj{x) = sin(27rx/L) + b, where b is now a 
dimensionless quantity. In our case, we have used b = 1.02. The behavior of the quantities of 
interest have been corroborated by Brownian dynamic simulations performed by integration 
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of the Langevin equation using the standard stochastic Euler-algorithm. The average particle 
current in x-direction has been derived from an ensemble-average of about 3- 10^ trajectories 
according to the following expression: 



In order to test the accuracy of the F-J description, we have evaluated the behavior of the 
nonlinear mobility as a function of the scaled force /, for different values of the periodicity L. 
In Fig. [2], the dependence of the scaled nonlinear mobility on / is depicted for two different 
periodic structures of spatial period L, and for h = 1.02. In both cases, the nonlinear mobility 
depicts a monotonic increase with increasing the force. The value of the scaling parameter 
/ = j3FL up to which the F-J approximation with spatial dependent diffusion coefficient 
D{x) provides an accurate description depends on the spatial period L. Consistent with 
the F-J approximation scheme, for rather smoothly varying cross-section, or equivalently 
large periods L, the agreement between our precise numerics and the analytic solution is 
attained for quite larger /-values as compared to the case with strongly winding structures 
or small periods. The analytic solution does not capture the correct limiting value for large 
/. Whereas the analytic result of the F-J equation tends to the value given by Eq. (fT6|) . the 
limiting value observed in the simulations always tends to 1 which corroborates with the 
deterministic limit. That indicates that in the case of very strong forcing, the particles almost 
travel ballistically, and thus do not feel the effects of the boundaries. This all is consistent 
with the accompanying breakdown of the F-J approximation in this strong forcing limit for 
which the deterministic dynamics dominates the transport. 

By comparing the exact numerical results with the analytic solution of the F-J equation, 
Eq. (IT^ . we can identify a critical value of the scaled force, /c, beyond which the relative 
error in the mobility exceeds a certain value. This critical force plays a similar role to that 
introduced in the previous section. 

For our example of a 2D channel whose shape is defined by Eq. ( 127|) . the validity criterion 
given by Eq. fl26l) simplifies to 
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FIG. 3: The dependence of the critical value of the scaling parameter, i.e. /c, on the periodicity 
L of a 2D channel defined by the dimensionless boundary function = sin(27rx/L) + 1.02 is 
depicted. For f < fc the relative error of the Fick-Jacobs description is less than 1%. The solid 
line is a L^-fit of the critical values obtained by comparison of the approximative analytic and the 
exact numerical results. 

thus predicting that the critical value of the force scales as L^. 

This prediction has been verified by the simulations. Fig. [3] shows the value of the critical 
force for a tolerance of 1% as a function of the periodicity L. The critical value of the force 
depends quadratically on the periodicity L^, as predicted. 

In Fig. m we illustrate, for the considered two-dimensional channel, the regions where 
the F-J approximation is accurate when compared with the simulations. We depict the 
dependence on the periodicity L of the maximal critical scaled force, obtained by comparing 
numerical results with the analytic solution for the particle current, for different required 
relative errors. This diagram shows the regions in parameter space in terms of / and L for 
which an accurate solution is provided. Thus, it is possible to provide an accurate result by 
using the analytic solution over a wide range of the scaled parameter and the periodicity. 



VI. THE HYPOTHESIS OF EQUILIBRATION IN THE TRANSVERSE DIREC- 
TION 

From our simulations, we can actually analyze the validity of the hypotheses of equilibra- 
tion in the transverse direction on which the F-J description relies. In Fig. [5l we show the 
steady state distribution of a certain number of Brownian particles for different values of 



12 



150 



100 



50 




10 12 



FIG. 4: The validity diagram of the Fick-Jacobs approximation for biased diffusion is obtained 
upon a comparison between the precise numerics with the approximative analytic solution, cf. 
Eq. (I14p for a 2D channel with boundary function uj{x) = sm{2iTx/L) + 1.02. The dependence of 
the critical value of the scaling parameter on the periodicity is depicted for three different relative 
errors; 1%: solid line, 5%: dashed line and 10%: dotted line. Below these limiting lines the analytic 
treatment agrees within the corresponding prescribed relative error. 



the scaled force /. As the force increases, we can clearly see how the particles are not homo- 
geneously distributed in y-direction, evidencing the failure of the equilibration assumption. 
This effect is specially dramatic in Fig. [5]^lower right panel), where the force is so strong 
that the particles do not fully explore the available space in y-direction. 

A more detailed analysis could be provided by checking the normalized steady-state 
probability distribution in the transverse direction at a given x-positions, i.e. 



(30) 



In Fig. [HI we represent the steady-state probability density at three different locations 
along the channel, corresponding to x/L = 0.2, x/L = 0.5, and x/L = 0.8. For small values 
of /, at x/L = 0.5, where lj'{x) = 0, the P^^iy) is flat, indicating a perfect equilibration in the 
transverse direction. However, at x/L = 0.2, and x/L = 0.8, the system is not equilibrated 
and P^^{y) is bell-shaped. Notice that for x/L = 0.2 and x/L = 0.8, the distributions are 
very much similar for small values of /, corroborating that the equilibration depends on 
C(j'(x)^, which is the same in both cases. 

At large force strengths particles concentrate along the axis of the channel. In this 
situation, the assumption of equilibration along the transverse direction fails. Particles would 
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FIG. 5: Steady-state particle distribution mapped into a single period of the 2D channel defined 
by the boundary function uj{x) = sin(27rx/L) + 1.02 with L = 1 for four different values of the 
dimensionless scaling parameter /: upper left panel: / = 0.2; upper right panel: / = 3.0; lower 
left panel: / = 7.0; lower right panel: / = 50. 
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FIG. 6: (color online) Normalized steady-state probability distribution of particles in y-direction 
PJ*(y), cf. Eq. (j30p for different values of the scaling parameter / - solid line (black): / = 0.2; 
dashed line (blue): / = 3.0; dotted line (green): / = 7.0; dashdotted line (red): / = 50.0, at various 
positions along the length of the channel with the boundary function lo{x) = sin(27rx/L) -|- 1.02, 
L = 1: x/L = 0.2 (panel (a)), x/L = 0.5 (panel (b)) and x/L = 0.8 (panel (c)). The grey regions 
indicate the outside of the channel. 



only feel the presence of the boundaries when they are close to the bottlenecks. Hence, in the 
limit of very large force values, the influence of the entropic barriers practically disappears. 
In this limit, the correction in the diffusion coefficient leading to a spatial dependency, i.e. 
D{x), overestimates the influence of the entropic barriers. 

Fig. [7] depicts the normalized stationary distribution function of particles in x-direction 
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FIG. 7: (color online) Normalized steady-state probability distribution of particles along the length 
of the channel P'^^{x\ cf. Eq. (|3ip . for the corresponding set of scaling parameters / as those chosen 
in Fig. [H 

at various force strengths, obtained from numerical simulations: 



P'\x) := /--^\\^' ""'^^ . (31) 



At very low force strengths particles are evenly distributed, and the probability distri- 
bution P^^{x) scales with the cross-section of the channel. When increasing the force, the 
maximum of P^^{x) is shifted towards the exit of the cell and the particles accumulate in 
front of the bottleneck. This behavior can also be observed in Fig. [51 In the large force 
regime P^^{x) is almost constant over a wide range of x- values, indicating a minor influence 
of the shape of the structure on the dynamics of the particles. In the limit of large /-values 
the plateau extends and covers the whole period. In this situation, a deterministic treatment 
of the problem leads to adequate results. 

VII. CONCLUSIONS 

With this work we have investigated the validity conditions under which the Fick- Jacobs 
approximation provides an accurate description of the biased diffusion of Brownian parti- 



cles in 2D and 3D confined systems. We have established a validity criterion 
terms of the sinuosity of the channel co'(x), as done in the unbiased case 
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brmulated in 
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and on the scaling parameter that causes forced diffusion. This scaling parameter com- 
pares the work done on a particle travelling a distance equal to the spatial period L of 
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the channel with the available thermal energy. Interestingly, the critical value of this scal- 
ing parameter up to which Fick- Jacobs equation holds depends on the square of the pe- 
riod. This dependence follows from the analysis of the different time scales that rule the 
biased, diffusive dynamics. We have constructed a validity diagram showing the region 
of parameters (spanned by / and L ) in which the Fick- Jacobs approximation describes 
the overdamped diffusive transport accurately cf. Fig. HI We have also investigated nu- 
merically the conditions for a fast equilibration in the transverse direction which is vital 
for the accurateness of the Fick- Jacobs approximation. The results presented evidence 
the usefulness of the Fick- Jacobs description with a spatially-dependent diffusion coeffi- 
cient at small applied bias. The obtainment of a simplified Fick- Jacobs like description 
that covers also for intermediate-to-strong values of the scaled force still presents an open 
challenge. The availability of such a description would be beneficial for the detailed under- 
standing of diffusive transport processes occurring in confined media, the latter dictating the 



nonequilibrium transport 
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3ehavior in a great variety of systems far from thermal equilibrium 
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APPENDIX A: PARTICLE CURRENT 



The evolution equation of the probability distribution of overdamped Brownian particles 
in ID can also be expressed as 



dt dx \ dx 
where D{x) is a position dependent diffusion coefficient, and A{x) is the free energy landscape 
for the reaction coordinate x. 
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Eq. (lAip results from the continuity equation: 



|p(x,t) = -^J(x,t), (A2) 



with the probabihty current: 



J(x, t)= - Dix) Q-PM-)^^m-) pu^ t) . (A3) 

ox 

In the case of a periodic diffusion coefficient, D{x + L) = D{x), and a tihed periodic free 
energy, A{x + L) = A{x) — FL, it is convenient to define the reduced probabihty density 
and the corresponding current. 



P(x,t) = ^P(nL + x,t), (A4) 

n 

J(x, t) = Y^ J{nL + X, t) , neZ. (A5) 

n 

By definition these functions are periodic with periodicity L, P{x + L,t) = P{x,t) and 
J{x + L,t) = J{x,t). Moreover P{x,t) and J{x,t) enter a continuity equation (1A2I) and 
P{x,t) is normahzed on any interval {x,x + L), provided that P{x,t) is normalized, e.g. 

dx P{x, t) = 1. In the steady state limit the probability current is a constant, J{x, t) — >■ 
J, thus Eq. ( 1A3I) becomes 

Multiplying both sides of Eq. (IA6p by 1/D(x)e^^*^'^\ and integrating over a period we 
obtain 

fx+L -| rx+L 



Jx D{x') J, dx' ^ > 



(AT) 



Using now the conditions P^^{x + L) = P^^{x) and A{x + L) = A{x)—PFL, Eq. (1A7I) reduces 
to 

J ^m-') dx' = P^\x) (1 - e-^^^) e'^^(^) . (A8) 

Jx D\X j 
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Rearranging the terms and integrating over to L, Eq. (lASP leads to 

J /^-'^^(^Mx r^''-i-e^^(^')dx' = (1 - e-^^^) r P'\x)dx, (A9) 
Jo Jx D{x') Jo 

which due to the normahzation yields 

J = ^ ^ . (AlO) 

/„^dxe-/^^W/r^dx'^e^^(^') 

The general relation between the stationary probability current and the steady state par- 
ticle current ((x)) is 



(x) = [ Jdx (All) 
Jo 

which implies (x) = J L. Then, the particle current reads: 

(x) = 5^ -J ^ (A12) 

\ / fL dx f3A{x) r^+^dx'^— e^^(^') 
Jo Jx Z)(x') 

Remarkably, this expression for the particle current is equivalent to the expression ob- 
tained via the mean- first-passage-time approach presented in Ref. [2^. 
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